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The Bethe-Salpeter equations for the quark-antiquark composite systems, qq, are written in terms 
of spectral integrals. For the qq-mesons characterized by the mass M, spin J, and radial quantum 
number n, the equations are presented for the following (n, M 2 )-trajectories: tyj, r/j, a.j, fj, pj, ujj, 
hj, and bj. 
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- - . I. INTRODUCTION 

o . 

The relativistic description of composite systems was always an actual and challenging task. The most frequently 
used technique, which takes into account the relativism of the constituents, is the Bethe-Salpeter equation [1]. But 
in the standard formulation of the Bethe-Salpeter equation, when the Feynman integration technique with mass-off- 
D ' shell amplitudes is used, one faces problems of the description of multiparticle channels and high-spin states. More 
appropriate technique for the high-spin composite systems is the dispersion relation approach, in particular, the most 

■ developed N/D- method [2]. However, our experience and intuition are based on the consideration of the potential type 
1 — I \ interactions, i.e., those which are associated with the particle-exchange mechanism. In terms of the 7V/D- method, 

one can easily relate the nearest left-hand side singularity of the V-function to the i-channel (or it-channel) meson- 
exchange diagram, but the reconstruction of the full set of left-hand singularities, when the interaction is given by the 
particle exchange or potential forces, is not a simple problem. Here we present the Bethe-Salpeter equation in terms 
of the spectral-integral technique which has advantages of both approaches discussed above: 
(i) in the spectral integrals, the mass-on-shell amplitudes are used, 
qq ■ (ii) the interaction terms are written in the potential or particle-exchange form. 

Moreover, in the spectral-integral technique one can use the energy-dependent forces as well. 

In [3-5], the dispersion- relation approach was applied to the description of the deuteron, being two-nucleon com- 
posite system, by treating nucleon-nucleon forces in terms of separable interactions. By using the interaction in a 

■ separable form one can work with mass-on-shell amplitudes and meson-exchange interactions. The expansion of the 
fS |' one-meson exchange diagrams in a series of separable interaction factors was developed in [6]. The principal points 

in the transformation of the standard Bethe-Salpeter equation to the dispersion-relation representation for the case 
of separable vertices were clarified in [7]. However, the representation of the meson-exchange diagram as a finite set 
of the separable vertices works successfully for mesons with nonzero mass only. For the long-range interaction, like 
confinement forces, the separable-vertex approach fails, thus forcing us to use not the standard V/D-method but the 
spectral-integral technique. 

The important ingredient of the dispersion relation and spectral-integration methods is the moment-operator ex- 
pansion. The elements of the moment-operator expansion technique were presented in [3,?, 8], and a systematic 
presentation of the technique may be found in [9]. 

The Bethe-Salpeter equation in the spectral integral representation is written here for quark-antiquark systems. 
Our attention is focused on the light-quark bound states, qq, where q = u, d, s:_for the simplicity's sake we consider 
here the systems built by quark and antiquark with equal masses: ud, uu, dd, ss. The treatment of heavy-quark 
composite systems, cc and bb, can be performed similarly. 

The necessity to deal with a full set of equations for the light-quark composite systems is governed by the rich 
information on the light-meson radiative decays that appeared recently [10-12] .The radiative decay data make it 
possible to restore the wave functions of mesons involved into reactions. Corresponding relativistic technique based 
on the consideration of the form-factor amplitudes in terms of the double spectral integrals was developed in [13] 
for pseudoscalar gg-mesons, and it was generalized for scalar and tensor qq-mesons in [14]. Finding out the meson 
wave functions in the spectral-integral form (or in the light-cone variables, see [13,?] for details) opens the way for 
direct reconstruction of the quark-antiquark interactions. The spectral-integral representation of the Bethe-Salpeter 
equation gives us the possibility to find out directly the interaction forces, provided the masses and wave functions 
of composite systems are known: this problem is discussed in Sect. 2 by using a simplified example of composite 
particles with spinless constituents. 

For the reconstruction of qq forces it is important for the light-quark qq-st&tes to lay on linear trajectories in the 
(n, M 2 ) plane, where n is the radial quantum number of the meson with mass M [15]. In more detail, the gg-statcs 
can be classified, within spectroscopic notations, as the n 2S+1 Lj levels, where S, L, and J refer to the spin, orbital, 



X 
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and total momenta, respectively. The analysis of spectra in the mass region 1950-2400 McV performed in [16] had 
fixed more than thirty new mesons which belong to the meson groups with positive and negative charge parities, 
(C = +/— ). Namely,for the (C = +)-states one has: 

1 Soqq — » 7r-mesons, 77-mesons, ry'-mesons, (1) 

1 D 2 qq — > 7r 2 -mesons, ?7 2 -mesons, 

3 Poqq — > a -mesons, / -mesons , 

3 P2qq — * a 2 -mesons, / 2 -mesons , 

3 Piqq — » a 1 -mesons, 

3 F 2 qq —> a 2 -mesons, / 2 -mesons , 

3 F 3 qq — > a 3 -mesons, 
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^4<7<7 - > a 4 -mesons , /4-mesons , 



and for the (C = — )-states: 



3 Siqq — > p-mesons, w-mesons , 0-mesons, (2) 
3 Diqq — > p-mcsons , 
3 D^qq — ► p 3 -mesons, 
P\qq —* /ii-mesons, 61-mesons , 
1 F^qq — > 63-mesons. 

The mesons measured in [16] as well as those accumulated in the compilation [17], being classified versus radial 
quantum number n, can be put, with sufficiently good accuracy, on linear trajectories in the (n, M 2 )-plane: 

M 2 = M 2 + /x 2 (n-l), n= 1,2,3,4,..., (3) 

with the universal slope /z 2 ~ 1.3 GeV 2 [15]. The linearity of trajectories, leading and daughter ones, was observed 
for the (J,M 2 ) plane too [15]. 

The linearity of trajectories on the (n,M 2 ) and (J,M 2 ) planes is in a good agreement with large-r behaviour of 
the confinement potential, V(r) ~ or, e.g., see [7], where the calculation of qq states from the groups (1) and (2) has 
been carried out. 

At the same time it is necessary to emphasize that for the low-mass states one can expect a violation of the trajectory 
linearity. For example, the 7r-meson is just an exception that is not surprising because of a particular role of the pion. 
The standard explanation is that the pion, being a low-mass particle, is determined by the instanton-induced forces, 
see [18,19] and references therein, although one cannot exclude an alternative modelling of the short-range forces. The 
problem of short-range forces is stressed by systematics of scalar states: the if-matrix analysis of tttt, KK, rjrj, and r]r( 
spectra [20] tells us that the lightest scalar-isoscalar state belongs to the flavor octet, but in the model calculations 
[18,19] the lightest state is close to the flavor singlet. We hope that a precise reconstruction of the qq forces can be 
facilitated by using the Bethe-Salpeter equation for the qq states in the spectral-integral form. 

So, we focus our attention on the reconstruction of the qq interaction, on the basis of the following triad: 

(1) the Bethe-Salpeter equation in the spectral-integral form, 

(2) the linearity of trajectories on the (n, M 2 ) and (J, M 2 ) planes, 

(3) the use of wave functions for low- and moderate-mass qq states found in the study of the meson radiative decays. 
The important point is that radiative decays can give us the information about meson wave functions which are 

now studied in the mass region 1000-1800 MeV: just the mesons from this region are determined by the short-r and 
intermediate-r forces, and only the forces from this r-region are not known sufficiently well thus being a subject of 
discussions and hypotheses. 

The paper is organized as follows. In Section 2 we recall basic statements of the Bethe-Salpeter equation written 
in terms of the standard Feynman diagram technique, give the elements of the dispersion relation iV/D-method, and 
clarify the interplay of these two methods by using separable vertices. The spectral integral representation of the 
Bethe-Salpeter equation is also written in this Section for the case of scalar constituents. 

In Sect. 3 the qq system is considered: the Bethe-Salpeter equations are written for the light-quark mesons which 
belong to the following (n, M 2 )-trajectories: 7rj, r]j, aj, fj, pj, loj, hj, and bj. 

In Appendices A, B, and C, the necessary auxiliary formulae arc presented which were used for deriving the 
equations. In Appendix D we collect equations which are rather cumbersome, these are the Bethe-Salpeter equations 
for oj, (f> 7 a 2 , and / 2 trajectories. 
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II. SCALAR CONSTITUENTS: DISPERSION RELATION 
METHOD AND THE BETHE-SALPETER EQUATION FOR COMPOSITE PARTICLES 



In this Section, we compare the Bethe-Salpeter equation for composite particles written with the use of Feynman 
diagrams with the equation in terms of the dispersion relations with separable vertices. This comparison gives us 
a guide for the transformation of Bethe-Salpeter equation with separable vertices into the spectral-integral Bethe- 
Salpcter equation with arbitrary meson-exchangc-type interaction. 

To simplify the consideration we deal here with scalar particles as constituents. 



A. Bethe— Salpeter equation in the Feynman-diagram technique 

Written in terms of the Feyman diagrams, the nonhomogeneous Bethe-Salpeter equation in the momentum repre- 
sentation reads: 

A(k 1 ,k 2 ;k';,k^ = v(k 1 ,k 2 ;k';,^)+ f Vikuhi-AX) 

J i(2irp 
S 4 (k[ +k' 2 -P) 



(m 2 -k[ 2 - i0)(m 2 - k' 2 - iO) 



A(k[,k 2 ;k'{k'>). (4) 



It is shown in Fig. 1 in the diagram form, and one can see there the notations for particle momenta. In (4) the 
constituents obey the momentum-conservation constraint: 



fci + k 2 = k[ + k' 2 = k'( + k' 2 ' = P, 

and V(ki,k 2 ;k[, k 2 ) is the irreducible kernel, i.e., the block without two-particle intermediate states (dashed block in 
Fig. 1). 

The scattering amplitude A(k\, k 2 ; k", k 2 ) determined by the Bethe-Salpeter equation (4) is the mass-off-shell 
amplitude. Even if we set k\ = k" 2 = k 2 = k 2 2 = m 2 in the left-hand side of Eq. (4), the right-hand side contains 
the amplitude A(k[, k' 2 ; k'{, k' 2 ') for k[ 2 ^ to 2 , k' 2 2 ^ m 2 . 

Let us draw the kernel V as a meson-exchange diagram; then, by iterating Eq. (4), we represent A(k\, k 2 ; k'{, k 2 ) 
as an infinite set of ladder diagrams of Fig. 2a. For further investigation it is important to fix intermediate states in 
the scattering amplitude. The ladder diagrams have two-particle intermediate states which can appear as real states 
at the cm. energies squared s = P 2 > 4TO 2 ,that corresponds to the cutting of ladder diagrams across the lines related 
to the constituents, see Fig. 2b. 

Such a two-particle state manifests itself as a singularity of the scattering amplitude at s = Am 2 . However, the 
amplitude A(k\, k 2 ; k'{ , k 2 ) considered as a function of s has not only this singularity but also an infinite set of 
singularities which correspond to the ladder-diagram cuts across meson lines associated with the forces: the example 
of such a cutting is shown in Fig. 2c. The diagrams, which appear after the cut procedure, are the meson-production 
diagrams, e.g., see Fig. 2d. 

So, in the complex s-plane the amplitude A(k\, k 2 ; k'{, k 2 ) has the following singularity: 

s = Am 2 , (5) 

which is related to the rescattering process. The other singularities are related to the meson production processes 
with cuts originating at 

s=(2m + nn) 2 , n= 1,2,3,.... (6) 
The four-point amplitude depends on six variables as follows: 

1,2 .2 >//2 h /l2 

s = (fci + k 2 f = (k'l + k'^f , (7) 
t= (fcj -k'l f = ik 2 -klf. 

The seventh variable, u = (k\ — k 2 ) 2 = (k'{ — k 2 ) 2 , is not independent because of the relation 

s + t + u = k 2 + k 2 + k'l 2 + kf. (8) 
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If the interaction creates a bound state, then the infinite set of ladder diagrams should produce the pole singularity 
in the amplitude. Near the pole, the scattering amplitude is determined by the diagrams of Fig. 3a type that means 
that in the graphical form the equation for composite system reads as Fig. 3b. In terms of the Feynman integral, it 
is as follows: 

A(k 1 ,k 2 ;P) = J ^y' 2 V(k u k 2 ;k[,k' 2 )x 

X (m 2 _ 4 ( -i0)(m2 P \q lQ) A ^ k '^ p y ( 9 ) 

The homogeneous Bethe-Salpeter equation (9), like nonhomogeneous one, works upon the mass-off-shell amplitudes; 
the multimeson production channels in (9) exist, and they are strongly related to the meson-exchange forces. 



B. Scattering amplitude in the dispersion relation N/D-method 

Let us summarize analytical properties of the discussed scattering amplitudes for two spinless particles (with the 
mass m) which interact through the exchange of another spinless particle (with the mass /i), Fig. 2a. This amplitude, 
A(s,t), has s- and i-channel singularities. In the t-planc, there are singularities at t — /i 2 , 4^i 2 , 9/i 2 , etc., which 
correspond to the one- or multi-meson exchanges. In the s-plane the amplitude has the singularity at s = Am 2 (elastic 
rescattering) and singularities at s = (2m + np,) 2 , where n = 1,2,... ,that corresponds to the production of n mesons 
with the mass fi. For the bound state with mass M, there exists a pole singularity at s = M 2 . If a mass of this bound 
state M > 2m, this is a resonance, and corresponding pole is located on the second sheet of the complex s-planc. 

The dispersion-relation N/D-method deals with partial-wave amplitudes. The s-channel partial amplitudes depend 
on s only. They have all the s-channcl singularities of A(s,t), namely, the right-hand-side singularities at s = M 2 , 
s = 4m 2 , s = (2m + /i) 2 , and so on, see Fig. 4. 

Left-hand-side singularities of the partial amplitudes are related to the i-channel singularities of A(s, t). The S*-wave 
partial amplitude is equal to 

l 

A{s)= J ^A(s,t(z)), (10) 
-l 

where t(z) — — 2(s/4 — m 2 )(l — z) and z = cos 9. The left-hand singularities correspond to 

t(z = -l) = (nn) 2 , (11) 

and they are located at s = Am 2 — fi 2 ,s = 4m 2 — 4^ 2 , and so on. 

The N/D-method provides us with the possibility to construct the relativistic two-particle scattering amplitude in 
the region of low and intermediate energies, where multiparticle production processes are not important; this region 
is shown in Fig. 4 by dashed line. If the threshold singularity at s = (2m, + /i) 2 is not strong (one-meson production 
process is suppressed), the region of partial amplitude under consideration can be expanded up to the next threshold. 

The unitarity condition for the partial-wave scattering amplitude (we consider the S- wave amplitude as an example) 
reads: 

Im A(s) = p(s) | A(s) | 2 . (12) 
Here p(s) is the two-particle phase space integrated at fixed s: 



p(s) = J d<S> 2 (P; k[, k' 2 ) = ±-^J S -^t, (13) 



d^ 2 (P;k' 1 ,k' 2 )= 1 -(2nrS\P-k' 1 -k' 2 )- ^ ""^ 



2 y J y 1 2 ' (2ir) 3 2k' 10 (2w) 3 2k' 2O 
In the TV/D-method the amplitude A(s) is represented as 
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where N(s) has left-hand singularities only, whereas D(s) has the right-hand ones only. So the TV-function is real in 
the physical region, s > 4m 2 . The unitarity condition can be rewritten as follows: 

Im D(s) = -p(s)N(s). (15) 

The solution of this equation is 

lX.)-i-/£«™.i-l*). 00 

4m 2 

In Eq. (16) we suppose that the CDD poles [21] are absent and we normalize TV(s) by the condition D(s) — > 1 as 
s — > oo. 

In principle, Eqs. (14), (16) provide us with a complete description of partial amplitude in the low-s region: the 
amplitude is determined by the A^-function being a set of the left-hand singularities which are due to the one-meson 
exchange (s = 4m 2 — fi 2 ), two-meson exchange (s = 4m 2 — 4^ 2 ), and so on. The right-hand singularities in Eqs. (14), 
(16) are uncoupled to the left-hand ones, opposite to the Feynman-diagram approach given by (4). It is important for 
the description of the realistic processes to have the left-hand and right-hand singularities uncoupled: a well-known 
example provides us with the pn-amplitude, with the deuteron quantum numbers, where the production of pions is 
suppressed (right-hand singularity at s = (2m + ^u^) 2 is weak), while the forces related to the pion exchange are 
significant (left-hand singularity at s = 4m 2 — fi 2 is strong) . 



C. N/D- method and separable interaction 

The N/D-method gives us the mass-on-shell partial amplitude, provided the TV-function is known. However, the 
TV-functions have rather intricate properties: they depend on the total number of the t- and u-channel exchanges and 
do not obey the factorization constraints, i.e., for different reactions the TV-functions may be independently different. 
As was stressed above, the spectral integral representation for the Bethe-Salpeter equation, keeping advantages of the 
dispersion relation method, is free from this problem: it uses t- and u- channel exchanges, with universal interaction 
blocks. 

As the first step in rewriting the Bcthc-Salpcter equation in the spectral-integral form, let us consider separable 
interaction as an example. For this purpose, we rewrite Eqs. (14), (16) introducing the vertex function 

g{s) = ^N{s). (17) 

Here we assume that N{s) is positive (the cases with negative TV(s) or with changing-sign N(s) need the introduction 
of several vertices). So, the partial wave amplitude A(s) written in terms of the separable vertex g(s) is given by the 
following series: 

A{a) = g(s) [1 + B(s) + B 2 (s) + B 3 (s) + • • •] g(s). (18) 

Its graphic interpretation is shown by Fig. 5. 

This set of diagrams can be rewritten in the form of the Bethe-Salpeter equation: 

oo 

A(s)=g 2 (s)+g(s) f ^g( s ')4^-A(s' , s), (19) 

J 7T 8—8 

im 2 

where A(s', s) is the energy-off-shell amplitude which enters the intermediate state of the diagrams of Fig. 5; one has 
for the energy-on-shcll amplitude A(s, s) — A(s).The interaction block is written as follows: 

V(s', S )=g( S ')g( S ) (20) 

Therefore, the Bethe-Salpeter equation is to be applied for the amplitude A(s,s'); it reads: 

oo 

A(s>, S )=g( S ')g( S )+g( S >) J ^ fl(a ")^^( s ", a). (21) 

4m 2 
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If the bound state exists, the amplitude contains a pole singularity at s = M 2 . Considering Eq. (21) near the pole 
and neglecting the non-pole terms, we have the following equation for the bound state vertex: 



oo 

G{ S \M 2 ) = g{s>) / ^ g ( s »)-£^ G (s",M% 



(22) 



4m 2 

where 

~A(s\ S )(M 2 -s) 



G{s',M 2 ) = 



G(s) 



(23) 

s^M 2 



The Bethe-Salpeter equation (22) gives us a guide for the consideration of general case, when the interaction is of 
the meson-exchange type. But before we need to consider in more detail the representation of the loop diagram. 

D. Loop diagram 

The loop diagram B(s) plays the decisive role for the two-meson amplitude, so let us compare in details the 
dispersion and Feynman integral expressions for B(s). 

Namely, the Feynman expression for Bp(s), with a special choice of separable interaction G(4fc^ + 4m 2 ), is proved 
to be equal to the dispersion integral representation, where the four-vector k± is defined as follows: 

2fc± = fci-fc 2 - ^~ 2 fc ' (fci + fc 2 ). (24) 

In this Section, we use the total-momentum vector P = k\ + k 2 , so it is convenient to write here P 2 but not s. 
The Feynman expression for the loop diagram reads: 



B (P 2 ) = J— f ( 4 ( Pfc 2) 2 AP 2 ~4k 2 +4m 2 ) 

F[ } (2n)HJ (m 2 - k 2 - iO)(m 2 - (P - k 2 ) 2 - iO)' { ' 

Since it is more convenient to treat composite system with light-cone variables, they are hereafter: 

k- = ^( fc 20 - k 2 z); k+ = -^=(k 20 + k 2z ); k 2 ± = k T ■ (26) 

We choose the reference frame, where Pt = 0. Then 

Pk 2 = P+k_ + P_k+, (27) 

and Eq. (25) takes the form 

- (in * < 28 > 

dk + dk-d 2 kr 



(2fc+fc_ -m% + iO)(P 2 - 2(P+fc_ + P_fc+) + 2fc+fc_ - m\ + iO) ' 

where m T = m 2 + k\ . It should be mentioned, that if G = 1, one can perform the integration over fc_ right now, 
closing the integration contour around the pole 

m% — iO 

k - - ~kr> ^ 

and we obtain the standard dispersion representation for the Feynman loop graph (x = k + /P + ): 
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1 



1 f 2 f {-2m) _ f ds 

(2^4 J dkT J dX 2x{P 2 - m 2 T /x - P 2 x + iO) ~ J w{s -P 2 - iO) X ( } 



o 



oo 

/dxdk\ 5{s — m\/[x{l — x)]) f ds ■ p{s) 
x{l-x) 16tt ~~ J ir{s - P 2 -i 



o)- 

4m 2 

The dispersion integral (30) is divergent at s — > oo due to G = 1, and it is the function G which makes Bp 
convergent in Eq. (25). Convergence of the integral (30) can be restored by the subtraction procedure. 

For G ^ 1, some additional steps are needed to obtain the dispersion representation, namely, we introduce new 
variables £ + and 

P+fc_ +P_fc+ = Vp 2 ^ 

P+fc_ - P_fc+ = VP2£-. (31) 
With these variables, Eq. (25) takes the form 

B F {P 2 ) = ' 



{2tt)H 

G 2 (4(£ 2 + m 2 T ))di + di-d 2 k T 



{£,+ -£,--m 2 T + ii)){P 2 - 2\AP2£ + + ^2 _ ^2 _ m 2 + iQ j 







= I dk 2 2d£_7rG 2 (4(£ 2 +m|)) x (32) 



(e + m 2 T ) + i0}[{£+ - VP 2 ) 2 - (e + m 2 ) + ,0)] 



The integration over £ + is performed by closing the integration contour in the upper half-plane, so the two poles, 
£ + = — y^; 2 . + TOy+zO and £ + = VP 2 — y^; 2 + m|,+i0, contribute. The introduction of a new variable s = 4(£ 2 +m|,) 
yields 



4m 2 

that is the dispersion representation of Eq. (16). 

So the hypothesis of the separable interaction gives us an opportunity to solve the Bethe-Salpeter equation easily. 
Within this hypothesis we can use different techniques: either Feynman integration, or spectral-integral representation, 
or light-cone variables. 



E. Spectral-integral representation and interaction forces 

The introduction of a separable interaction is not the only way to make the Bethe-Salpeter equation easily solvable. 
The main point in handling the Bethe-Salpeter equation is to control the right-hand side singularities, especially those 
related to multimeson production, at s = {2m + p) 2 , s = (2m + 2/i) 2 ,... , and it is the spectral integration technique 
which enables us to control the multimeson production processes. 

The spectral integral representation is based on the following corner stones: 
(i) constituent particles in the intermediate states are mass-on-shell {k 12 = m 2 and k' 2 = m 2 in Fig. 1); (ii) there is 
no energy conservation in the interaction processes (s s' s" in Fig. 5). 

Based on these statements we consider potential interaction, or the particle-exchange interaction, by using the 
spectral-integral diagrams. Consider as an example the interaction associated with the t-channel exchange by a 
meson with the mass p: 

V{k 1 ,k[) = ^— t , t={k 1 -k[) 2 . (34) 
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In the cm. system, which is the most convenient for the consideration, one has for the four-momenta of the 
constituent particles: 



fa = (k ,k) = hl±,nJ S --w?y (35) 



^1 — (^0' & ) — 

where n 2 = n' 2 = 1. The interaction block V(ki, k[) can be expanded in a series with respect to z — (nn 1 ). In this 
way, we may obtain the interaction for different partial waves. For example, the interaction in the wave with the 
angular momentum L — is equal to: 

l 

Vo(s,s') = J yV(fcl.fci). (36) 

-1 

Actually, Eqs. (34), (35), (36) allow us to generalize the procedure with the separable interaction considered above. 
Indeed, expanding (36) in a series with respect to orthogonal functions, one has 

Vo(*,*') = $>n(*)0n(*')> (37) 

n 

that is a separable interaction in a generalized form, supposing the choice of the functions allows one to use a 
restricted number of terms in (37). Separable interaction taken in such a form was used in [3,4] for the description of 
nucleon-nucleon interactions by considering the deuteron within dispersion-relation technique. 




F. Spectral-integral representation of the Bethe— 
Salpeter equation for composite system 

First, we consider the case of L = for scalar constituents with equal masses, though not identical. The bound 
system is treated as a composite system of these constituents. Furthermore the case L ^ is considered in detail. 



1. Bethe- Salpeter equation for vertex function with L — 
The equation for the vertex composite system — > constituents, shown graphically in Fig. 3b, reads: 

oo 

G ( s ) = / f d$ 2 {P'-X,k' 2 )V{k u k 2 -XX) g , ^M*- iQ ' m 



4m 2 



where the phase space is determined by Eq. (13). Scalar constituents are supposed to be not identical, so we do not 
write additional identity factor 1/2 in the phase space. 

The equation (38) written in the spectral-representation form deals with the energy off-shell states s' — (k' 1 + k / 2 ) 2 ^ 
M 2 , s = (ki+k 2 ) 2 ^ M 2 and s ^ s'; the constituents are mass-on-shell, k' 2 = m 2 and k' 2 2 — m 2 . We can use alternative 
expression for the phase space: 

rf$ 2 (P'; k[,k' 2 ) = p( S ')y = d$(k>), z = ~^=, (39) 

where k = (fa - fc 2 )/2 and kl = (k[ - k' 2 )/2. Then 

OO 

G(s) = / V / M{k ' ] V (S ' S '' {kk ' )] s> G M 2 - i - (40) 

4m 2 

In the cm. system (kk 1 ) = -(kk') and Vk 2 = V^k 2 = i|k| and Vk 12 = V^k 72 = »|k'| so z = (kk')/(|k||k'|). The 
phase space and spectral integrations can be written as follows: 
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oo 

f ds' f , ,,,,,,, f dk' 

Am 2 



(41) 



where k' = V to 2 + k' 2 . In the cm. system Eq. (38) reads 



,2. Bethe-Salpeter equation for the (L — 0)-wave function 
Now consider the wave function of composite system, 



To this aim, the identity transformation upon the equation (42) should be done as follows: 

oo 

(s - M2) I^& = / T / d<J>(fc,)t/ (s ' s ' (fcfc,)) 5^ ' (44) 

4m 2 

Using the wave functions, the equation (44) can be written as follows: 

oo 

{s-M 2 )^{s) = J ^- J d$ 2 {k')V (s,s' ,{kk'))^{s') . (45) 

4m 2 

Finally, using k' 2 and k 2 instead of s' and s, 

1>(8) -> V(k 2 ), 

we have: 

dk' 



^-^U( S , S ',-(kk'))^(k' 2 ). 



(46) 



This is a basic equation for the set of states with L = 0. The set is formed by the levels with different radial excitations 
n=l,2,3,..., and relevant wave functions arc as follows: 

?Ai(k 2 ), ^ 2 (k 2 ), ^ 3 (k 2 ),... 

The wave functions are normalized and orthogonal to each other. The normalization/orthogonality condition reads: 

/r/k 
^-^-^(k 2 )^(k 2 ) = ^. (47) 

Here S nn ' is the Kronecker symbol. The equation (47) is due to the consideration of the charge form factors of 
composite systems with the gauge- invariance requirement imposed, see for detail [13]. This normalization/ 
orthogonality condition looks as in quantum mechanics. 

Therefore, the Bethe-Salpeter equation for the S'-wave mesons reads: 



oo 

.'2 



r die 12 

4(k 2 + m 2 )^(k 2 )-y ^U o (k 2 ,k' 2 )0(k' 2 )^(k' 2 )=M 2 V„(k 2 ) , (48) 



where 
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*< k ' 2 > = rM ■ «•) 

The f/'nlk 2 ) presents a full set of wave functions which are orthogonal and normalized: 

oo 

r dv 2 

J — ^(k 2 )0(k 2 )^(k 2 ) = <5 ab . (50) 
o 

The function y (k 2 ,k' 2 ) is the projection of potential V(s,s', (kk 1 )) on the S*-wave: 

V (^) = J d ^J d ^V(s,s>,-(W)). (51) 

Let us expand Vo(k 2 ,k' 2 ) with respect to full set of wave functions: 

yo(k 2 ,k' 2 )=^^(k 2 ),i° b V b (k' 2 ), (52) 

a, b 

where numerical coefficients are defined by the inverse transformation as follows: 

oo 

v (0) = J ^!^! V , a(k 2 )0(k 2 )yo(k 2 )k/ 2 )0(k/2)v , b(k ,2 ) _ (53) 


Taking account of a series (52), the equation (48) is rewritten as follows: 

4(k 2 + m 2 )^„(k 2 ) - ^(k>W = M 2 ^„(k 2 ). (54) 

a 

Such a transformation should be carried out upon the kinetic-energy term, it is also expanded in a series with respect 
to a full set of wave functions: 

4(k 2 + m 2 )V.„(k 2 ) = ]T i^a(k 2 ) , (55) 

a 

where 

oo 

K na = J _^ Q (k 2 )0(k 2 ) 4(k 2 + m 2 )^„(k 2 ) . (56) 
o 

Finally, the Bethe-Salpeter equation takes the form: 

]T i^ a (k 2 ) - Y, «iaVa(k 2 ) = M 2 ^„(k 2 ) . (57) 

a a 

We take into account that Vnl = i>cto- 

The equation (57) is a standard homogeneous equation: 

^ Sna ^ a (k 2 ) = M 2 ^„(k 2 ) , (58) 

a 

with s na = K na — Vnd ■ The values M 2 are defined as zeros of the determinant 

det\s- M 2 I\ = , (59) 

where / is the unity matrix. 
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3. The Bethe-Salpeter equation for the states with arbitary angular momentum L 

,,(£) 



For the wave function with arbitrary angular momentum W„) Ml p t ( s )' we use * ne following ansatz: 

b (L) 

{n)ni,...,iXL 



<L,... )W («) = 4t ) ..^(^ ) w- (so) 



The momentum operator xj£l.. tttL (k) was introduced in [9], we remind its features in Appendix A. 
The Bethe-Salpeter equation for the (L, n)-state, presented in the form similar to (48), reads: 



x 



4(k 2 + m 2 )AW..^(fc)^)(k 2 ) - X%l.,» L {k)x 

OO 

J ^VUs,s')Xl(k'^(^)^^ 2 ) = M 2 AW..,, L (fc)^ L) (k 2 ), (61) 



where 



xl{k , 2) = Jd^_ (xW ^y = a(L)(k'*) L = a{L)(-^)\ (62) 

a{L) = (2L_m a(0) = 1 (63) 

The potential is expanded in a series with respect to the product of operators xj^... tllL (k)xjff... tflL (k'), that is, 

V(s,s',(kk'))= J2 ^t.,^)VL(s,s')X^ L (k'), 

L,ni_...n L 

Xl{k 2 )V L { S ,s')Xl(k' 2 )= (64) 

Therefore, formula (61) reads as follows: 

OO 

^Vi( S , S ')«W(-k' 2 ) L ^(k' 2 )^)(k' 2 ) = (65) 

o 

As compared to (48) this equation contains additional factor Xl(k' 2 ); still, the same factor is in the normalization 
condition, so it would be reasonable to insert it into the phase space. Finally, we have: 

OO 

4(k 2 +m 2 )V^(k 2 )-/ ^^( S , S ')^(k' 2 )^(k' 2 )^M 2 ^)(k 2 ) , (66) 
o 

where 

L (k' 2 )- a (L)(k' 2 ) L 0(k' 2 ), V L ( S , S ') = (-l) L V L (s, S '). (67) 
The normalization condition for a set of wave functions with orbital momentum L reads: 

OO 

J ^^)(k 2 )^(k 2 )# ) (k 2 )= < 5 Qb . (68) 
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One can see that it is similar to the case of L = 0, the only difference consists in the redefinition of the phase space 
4> — > 4>l. The Bethe-Salpeter equation reads: 



with 



^e o )e)(k 2 ) = M^i L) (k 2 ), 



(69) 



(L) =R (L) _ (L) 

oo 



oo 

^e ) (k 2 )^(k)4(k 2 +m 2 )e)(k 2 ). 



(70) 



Using radial-excitation levels one can reconstruct the potential in the L-wave, then to reconstruct, with the help of 
(64), the i-dependent potential. 



III. QUARK-ANTIQUARK COMPOSITE SYSTEMS 

For the qq system, the Bethe-Salpeter equation for the wave function with the total momentum J, angular mo- 
mentum L = | J — S\ and quark-antiquark spin S can be conventionally written as follows: 



(*-^)*!3^*>=y^ 



d 3 k' 



-V(s,s',(kk'))§ 



(S.L.J) 



(*') , 



where 



k=^(h-k 2 ), 3 = {k l +k 2 ) 2 , k' = \{k[-k' 2 ), s' = (k[ + k' 2 f. 



The wave-function operator with the fixed quantum numbers is presented as 

where Q is the moment operator for the qq system. 
The potential operator can be decomposed as follows: 

V (s, s', {kk')) = V} 0) («, a'. (**')) 0/ ® Oj , 



(71) 

(72) 
(73) 

(74) 



where I — S,V, T, A, P is a full set of Dirac matrices in the t-channel: 

Oi = I, 7 M , icr^, »7m75, 75- (75) 
The potential operator V (s, s', (kk')) can be decomposed in the s channel by using the Fierz transformation: 

V (s, s', (kk')) =Y,J2 ^/ (0) ( s > s '> ( kk ')) °ic (Oc ® Oc), (76) 



7 c 



where Cj c are coefficients of the Fierz matrix: 



C 



Ic 



(\ 


i 

4 


1 

8 


1 1 \ 

4 4 ' 


1 


1 
2 





1 "I 


3 





1 

2 


3 


1 


1 

2 
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(77) 
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Denoting 



we have 



V c (a, a', (kk 1 )) = {a, a', (kk 1 )) C Ic , 

i 



V {a, a\ {kk')) = J2(O c ® O c ) V c {a, a\ {kk')) = 



= (I <8> I) V s (a, a', (kk')) + O 7„) W (s, (fcfc')) + (tCT„„ <8> io^ v ) x 

x V T (s, a', {kk')) + (i7„75 ® ^75 ) Va (s, (&&')) + (75 <8> 7s) Vp (s, a', {kk')) . 

Let us multiply Eq. (71) by the operator Q < f 1 [ L .' l f) (k) and convolute over the spin-momentum indeces. 
definition V(s, a', (kk')) -> (fcj + m')V(a, a', (kk'))(-k' 2 + m') one has: 

{a M 2 ) Sp pg^/fcXfci + m)Q^)(k)(-k 2 + m) 
= J2Sp[0 c (k 1 + m )Qit L :, , }(k)(-k 2 +m)\ J ^-^V c (s, 3 ',(kk'))x 
xSp [(%[ + m')O c (-k' 2 + m'^f^J^k') . 

Here we define to' as k'\ = k'\ = to' 2 . We have four states with the qq spins 5 = and S = 1: 

1) S = 0; L = J, 

2) S = 1; L = J+ 1, J, J - 1. 

These states are constructed from the operators [9] as follows: 

Q { °dil(k) = i 75 x^ j (k), 
QiY^; J \k) = ^x^l a (k), 

Q < ili'. J .fj, 1 / J Hk) = 7„ ^j ^k). 

For these operators, the wave functions read: 
for J = L + 1, 



for J = L — 1, 



for 5= 1, L = J± 1, J, 



*!^^( fc ) = ^' J) w4 1 ' J - 1 - J) (fc 2 ), 



^(l (7±l),J) (fc) = , ^xj-x.j, ( fc ) + R.tfJY™ (*). 



(i,j+i,j) 



where Aj and are the mixing coefficients with j = 1,2. 
These wave functions arc normalized as follows: 



/ 



rf 3 fc 
(27r) 3 fc 



(-l)Sp 



^(S",Z/, ,J') 
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A. Equation for (5 = 0, J = i)-state 

The equation for the state with 5 = 0, J = L reads: 

(a - M 2 ) xW. w (fe) Sp \i lb (ki + m)i l5 (-k 2 + m) 



(89) 



xXtf..^ (k) ^ JJ) (k 2 ) = Xj/l^k) SP [f c (fci + m)i l5 (-k 2 + m) 



I 



d 3 k' 



T V C (s, s', (kk')) Sp i l5 (k[ + m')F c (-k' 2 + m') 



Now consider the left-hand side of the Eq. (89). Using the traces presented in Appendix B and convolution of 
operators from Appendix C, we have: 

XlH^(k)Sp [nsCfci + m)i l5 (-k 2 + m)l X^(k) = 



= X^ik) (-2 S )4f).. w (fc) = -2sa(J) k 2J . 
The right-hand side of the equation is calculated in two steps: first, we summarize with respect to c: 

A(s,s',(kk'))= A c ( s ,s',(kk'))V c ( Sl s',(kk')) = 



c=T.A.P 



(90) 



Y Sp Fc (h + m)i-y 5 (-k 2 + m) Sp ij 5 (k[ + m')F c (-k' 2 + m! 



c=T,A,P 



xV c (s,s',(kk')). 

In Appendix B the trace calculations are presented, and the values 
A c (s, s', (kk')) are given. In this way, the sum is written as follows: 

A{s,s',{kk'))= A c (s,s',(kk'))V c (s,s',(kk')) 

c=T,A,P 



— 4vsa' X 



(91) 



G^Vp (s, a', {kk')) + 4mm' V A (a, s', (kk')) + 8(kk') V T (a, s', (kk')) 

At the second step, the convolution of operators is performed by using equations of Appendix C and recurrent formulae 
for the Legendre polynomials: 

that allows us to write the Bcthe-Salpetcr equation in terms of the Legendre polynomials (recall that z = 
(kk') / (Vti^Vk^ 2 ) and Vk 2 = iy/s/A — to 2 , Vk^ = i^Js' '/4 — to' 2 ). For the exceptional case J = we put P-i(z) = 
As a result we get: 



4f Uj ( k ) A ^ i kk ')) X l J l^ ( fc ') = «( J ) (v^V^) 3 (-4^) x 



(92) 
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8^±^V¥Vk^Pj+i(z) Vt (a, a', {kk')) + 



+ (y7^Vp(s, S ',(kk')) + Am m 'V A ( Sl s',(kk'))^ Pj(z) + 
J 



-8- 



2J+ 1 

Substituting the obtained expressions into Eq. (89), we obtain: 



fc 2 v7c' 2 P 7 _i(z) V T (a, {kk')) 



(s-M 2 ) (-2s)a(J)k 2J ^ J \k 2 )= j 



d 3 k' 



(2nfk' 



(93) 



x (-4V^F) a( J) (VPVfc 72 ) ' 



*2J+ 1 



k 2 Vk' 2 P J+1 (z) V T (a, a', (kk')) + 



+ (y/rfVp (a, a', (kk')) + 4mm' V A (a, a', (fcfc'))) Pj(z)+ 



J 



+8 2J+T V^Vk'ZPj-tiz) V T (a, a', (kk')) 
Expanding the interaction block in the Legendre polynomial series, 

V c (s,s',(kk')) = Y / Vy ) (a,s')Pj(z) = 



¥°' J ' J) (k' 2 ). 



(94) 



= V c J) ( s ' s ') a ( J ) (-vWjfc*) J Pj(z), 
j 

and integrating over angle variables in the right-hand side by taking into account the standard normalization condition 
f\ dz/2 P 2 j(z) = 1/(2 J + 1), we have finally: 



OO . 

(a-M 2 )^\a) = j ^ P ( s ')(-k' 2 y2f- 



(95) 



J+l 



4m' 2 



" 8 2TTT €(J + 1)fc2fc ' 2 V ? J+1) (s ' s ' } + Vss ' €(J) ^ (s ' s ' } 



+4mm' £(J) (s, S ') - 8^^ £(J - l)^ 7 " 1 ) («, *') 
£ y Pj(z)K («, (kk')) = £ ^Pj(z) V C (J,) (*, s') Pr(z)a(J') x 

x (- Vk^j J = v} J) (s, s') (- VP v^) J = 



where: 



£(J)K (J) (*,*') (-V^V^)' 

a(J) (2J-1)!! 



2J+ 1 (2J+ 1) • J! ' 



(96) 
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1. Equation for the pion (M 2 ,n) -trajectory 



The pion states which belong to the pion (M 2 , n)-trajectory obey the following equation: 



ds' 



(97) 



4m 2 



-\k 2 k' 2 (s, s') + V^V { p ] (s, s') + Am 2 vf ] (s, s') 



Recall that n is the radial quantum number, and the following states with different n are located on the discussed 
(M 2 ,n)-trajectory: tt(140) with n = 1, tt(1300) with n = 2, tt(1800) with n = 3, tt(2070) with n = 4, tt(2360) with 
n = 5, and so on. 

The wave functions of the states laying on the (M 2 , n)-trajectory satisfy the orthogonality/normalization constraint: 



oo 



(98) 



4m 2 



The factor 2s is due to summing over the spin variables of quarks. 

Expanding the interaction block over full set of the radial wave functions, we can transform (97) into a system of 
the linear equations of the type of (58). 

2. Equation for the r\ (M 2 ,n) -trajectory 

The 77-states have two components, nn = (uu + dd)/y/2 and ss. We write rj n = cos9„ nn + sinO„ ss, and 
n' n = — sin 6„ nh + cos 0„ ss. For the lightest mesons 7/(550) and r/(958), one has cos 9i =~ 0.8 and sin 61 ~ —0.6. 
Correspondingly, we have two equations for the wave functions which describe the nn and ss components: 



(-M 2 )C»( S )-e n 



- P{ s') 2^- x 



(99) 



4 m 2 



.°k 2 k ,2 v w 



(s , S ') + ^V^ nnhP (s, s') + *m 2 V^ nn)tA (s, s') 



x<(S„( s ')cose n + / ^ Ps (s') 



ds' 



4m? 



3^ ft a v (sg— nft),T 



(*,*') + 



(o) 



',(0,0,0) , m 



where p s (s') refers to the ss phase space and y/kf 2 — i\J s' /4 — m 2 . The second equation, for the ss-componcnt, 
reads: 



(^M 2 )<ft 0) „(s)sine„^ J ^ S ') 2 V7 



(100) 



4m 2 



--/r 2 A-' 2 ^ (1) 



. sS)jT (s,s') + Vs^<1_ 



► ss),P 



(s, s') + 4mm s V, 



(o) 

(nn— >ss),A 



(s,s') 



16 



ds' 



{8,8') 



4 m. 2 



+^<L ss -),p + 4m? V^ sShA (8, s')] ^g n (8') sin 6, 



The wave functions ib 1 " ,' ! ? (s) and tb {u /-l'' (s) satisfy the normalization condition within an obvious change of the 

integration region for the ss component: 4m 2 — ► 4m 2 . 

The following states are located on the -q and rf (M 2 , n)-trajectories: 

1) ry-trajectory: ?7(550) with n = 1, ry(1295) with n = 2, r/(1700) with n = 3, r/(2010) with n = 4, 7,(2320) with n = 5, 
and so on. 

2) //-trajectory: 7/(958) with n = 1, ?7(1440) with n = 2, 7/(1820) with n = 3, and so on. 



',(0,0,0) 



3. Equation for the 61 (M 2 ,n)-trajectory 
The equation for the states with 5 = 0, L = l, J=l reads: 



(-M 2 )<:^( s) = -/^( S ')fc' 2 ^ 

4m 2 



(101) 



24 



y fc 2 fc' 2 t4 2) (s, s') + \/s7 v£ 1} («,«') + 4 ™ 2 ^ ( s > s ') - 8 v t } ( s > s ') 



The following states are located on the b\ (M 2 , n)-trajectories: &i(1235) with n = 1, &i(1640) at n = 2, &i(1970) with 
n = 3, 6i(2210) with n = 4, and so on. 

The wave functions of the &i-states laying on the (M 2 , n)-trajectory satisfy the orthogonality/normalization con- 
straint: 



oc 

/ %{s)28k 2 ^{8)^\s)=5 nl , n . 

4m 2 



(102) 



4- Equation for the hi (M 2 ,n)-trajectory 

The /ii-states have two components, nh — (uu + dd)/>/2 and ss; we write /ii,„ = cos9„ nn + sin9„ ss. Corre- 
spondingly, we have two equations for the wave- function nn-component: 



4m 2 



(103) 



-^fc 2 fc' 2 y/ 2) 
5 



_s „ (s, s') + \fss 1 v!' 1 l -, p (s,s') + 

(nn^nn),T \ ' ' (nn—>nn),P \ ' ' 1 



+4m 2 V, (1 > , (s, a') - 8y, W (s, a')l V'w'H («') cos9„- 

(nn-tnn),A V ' / (nn—>nn),T \ ' / ' n 
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-/ 

4m2 



ds , , io 2 is' 



24 



n, n. s ( ss _ mn ) ) T' 



(S,S') + VSS>V ( V, P {8,8') + 



(ss^nn),P 



+ Amm s V$_. nfiM (s, s') 8 V^ nfi)iT (s, s')} sin Q n , 



where p s (s') refers to the ss phase space and yk^ = zyV/4 — m^. For the ss-component we have: 

oo 



4m 2 



(104) 



24, 



_±ifc 2 fc' 2 v/ 2 l T ( s , s') + V^v ( [L l p (s, s') + 



(i) 



5 s (nrt— »ss),T 



+4mm s ^L-.-),^ ( s > s ') - 8 ^U,),t («> «') 



(0) 



oc 

-/ 

4m 2 



ds . ,o 2 / s' 
7T 3 V s 



24 



+4m 2 , (s, a') - 8 V, w T (a, a' 



(0) 



^M 1 ; s 1 ) ) ,n( s ')sine„, 



The wave functions V^'(nn) n( s ) anc ^ ^^'(s's) n( s ) sa tisfy the normalization condition within the obvious change the 
integration region for the ss component: 4m 2 — > 4m 2 . 

The following states are located on the hi (M 2 , n)-trajectories: 

1) /ii(1170) with n = 1, h^WOO) with n = 2, ^(2000) with n = 3, hi (2270) with n = 4, and so on; 

2) hi (1390) with n = 1, hi (1780) with n = 2, hi (2120) with n = 3, and so on. 



5. Equation for the ir 2 (M 2 ,n) -trajectories 
The equation for the ^-states (5 = 0, L = 2, J = 2) reads: 

oo . 

(s-M 2 )^(s) = / ^ P ( 3 >)k>* 3 -f-x 

4m 2 



(105) 



-yfc 2 fc' 2 V^ 3) (a, s') + ^V { P 2) (a, s') + 4m 2 Fj 2) (s, s') - ^V™ (a, s') 



The following states are located on the ir 2 (M 2 , n)-trajectory [15]: 7r 2 (1670) with n = 1, 7r 2 (2005) with n = 2, 7r 2 (2245) 
with n = 3, and so on. 

The wave functions of the 7r 2 -states satisfy the orthogonality/normalization constraint: 



oo 

/ ^P(*)fc 4 2 S a(2)^°; 2 ; 2) ( s )<„ 2 ^( s ) = ^,„. 

4m 2 



(106) 



where a(2) is determined by Eq. (63). 
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6. Equation for the 772 (M 2 ,n) -trajectory 



The ^-states have two components, nn — (uu + dd)/y/2 and ss. We write 772, n = cos8„ fin + sin6„ ss, and, 
correspondingly, we have two equations for the wave functions: 



00 . 

(--^£$U->«» e » = / v^' )fc 'iv7 

4m 2 



(107) 



40 , 



^k* V<2> T (s, s>) + V^V%> p (s, s>) 



(nn—>nn),T 



(2) 

(nn—>-nn),P 



+4™ 2 V^ nn)tA (s, s') - f l^nn )fT (a, 



00 
4ra 2 



3 /£ 
5 V s 



*j s (ss~*nn),T 



(s,s') + 



+4mm s V, 



(2) 



„^ _ 32 t7(i) 



V'i°'S ) ,„( s ')sine„, 



and 



ds' 



3 /s 



4m 2 



(108) 



40 



(nn— >ss),P 



(s,s') + 



+Amm s V { ^ sShA ( S , S ')-^V ( il) 



g (nn^ss),T 



(s,s') 



4m 2 



g (ss-^ss),T 



{8,8') 



('0 2 2') fO 2 2) 

The wave functions W^r^n) «( s ) anc * '(«'») n( s ) sa tisfy the normalization condition (106), with obvious change of the 
integration region for the ss component: 4m 2 — > 4m 2 . 

The following states are located on the 772 (M 2 , n)-trajectory [15]: 77 2 (1645) with n = 1, 77 2 (2030) with n = 2, 
?72(2250) with n = 3, and so on. 



19 



7. Equation for the bz (M 2 ,n) '-trajectory 



The equation for the 63-mesons (S = 0, L = 3, J = 3) is as follows: 

r2\ ,/,(0,3,3), o x _ 



(-M 2 )<;f)( s ) = -/^P( S ')^^ 

4m 2 



(109) 



--^fc 2 fc' 2 V4 2) ( S , s') + V^V^ («, s') + im'V^ (s, s') gv4 0) 



According to [15], the following states are located on the 63 (M 2 , n)-trajectory in the mass region below 2400 MeV: 
& 3 (2020) with n = 1, 6 3 (2245) with n = 2, and so on. 

The wave functions of the states laying on the (M 2 , n)-trajectory satisfy the orthogonality/normalization constraint: 



4m 2 



The factor a(3) is given by the Eq. (63). 



(110) 



B. Equation for the (5 = 1, J = L)-state 

The equation for the (5 = 1, J = L) state reads: 

(s - M 2 ) ep Vl w 3 P Vl Z < £l 1 ... flJtV3 {k)Sp 7^(fci + m)j^(-%2 + m) 



(111) 



x^ 1 ^s^ 1 4 J 2 1 - w ,«s( fc )^ 1 " 7,J) ( a ) = 

= ep> VlV2V3 P Vl zii] ll ... llJtV3 {k) S P [Fc (ki + m)7^(-fc 2 + m) 



(2tt)^ 



7 K {a, a', (kk')) Sp ri(k[ + m')F c (-k' 2 + m') 



The left-hand side of the equation is calculated by using the trace and operator convolutions given in Appendices B 
and C: 



e/ 3vi v 2 v 3 Pv 1 Zl 2 j ll ... llJ , V3 (k)Sp |r a (h +m)^(-k 2 +m) 
^^P^l.^^k) = -2s 2 + f a(J)k 2 -'. 



(J + 1)2 



As before, the right-hand side is calculated in two steps. 
1) We calculate traces: 



Bp, a ,{s,8 l ,{kk'))= J2 {Bc)f3> a >{s,s',{kk'))V c {s,s',{kk')) 
c=T,A,V,S 



(112) 



(113) 
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= J2 S P & + m )^' (-^ + m )] 5 P (*i + m ')^c {-% + m') 

c=T,A,V,S 

xV c (s,s',(kk')). 
Following the items presented in Appendix B, we write: 

Bp a , (a, a', {kk')) = g^, a , 4VsV \V^V V (s, a', {kk')) + 



(114) 



+8mm' V T (a, a', {kk')) + Wk^Vk' 2 z V A (s, a', {kk')) + 
+QAmm'kj,k'£ V s (a, a', (fcfc')) - 16kfik£, V7?V A (s, s', {kk')) + 



-16 



s'fc^fci + sk'^k'Jr + 4 z Vk*Vk*kfrk£ 



V v (s,s',{kk')). 



2) The convolutions of the trace factor -Ba< Q ' (a, a', {kk')) with angular momentum wave functions are presented in 
Appendix C; we have: 



E-fi r v\V2V3 Zv 2 ^ ll ■■■^ij {k) Bf3' a ' (a, S , (/c/l )) X 

J . A ,.J(2J + 3) 2 



xe a ^ 3 P il Z^l... tijiS3 {k') = a{J) (y&Vw*) {-Ass')- 



(J+l)3 



+ 



j+i 



+8mmV T (a, a', {kk 1 ))) Pj{z) + A——sJk 2 ^k' 2 Pj^{z)V A (a, s', {kk')) 
Inserting these expressions into Eq. (Ill), we obtain: 

(a-M 2 ) (-2a 2 )fc"v4 W) (*) = / ^(-W)(^Vf)''x 

t-jL-V^V^Pj+i^Vk (a, s', {kk')) + sf^V v (a, s', {kk')) Pj{z)+ 

+8mm'V T (a, a', {kk')) Pj{z) + A^^Vk^Vk^Pj^{z)V A (a, a', {kk')) 

x^' J ' J) {s'). 

Expanding the interaction block according to (94) and integrating both sides over dz/2, we get: 



oo 

(a-M 2 ) ^ J \s)= I ^p{ s '){- k ' 2 y2 S - 

J 7T S 



J 



4m' 2 



4 27TT aJ + 1) k2k ' 2 v * +1) (s ' s ' } + Vss ' €( J) v ^ ' (s ' s ' )+ 



+8mm' £( J) V± J) (a, a') - 4^±1 £( J - 1) v]- 7 " 1 ) (a, ^ 
Normalization condition for the {S = 1,J = L) wave functions reads: 



^> J ' J) (s>). 



OO 

J d i p{s)k2J 2s2,I jjTW a{J) ^ JJ) ^ JJ) ^ - <w- 

4m 2 



(115) 



(116) 



(117) 



(118) 
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1. Equation for the ai (M 2 ,n) -trajectory 
The ai states (S = 1, L = 1, J = 1) obey the Bethe-Salpeter equation: 



( S -M 2 )^)( S ) 



(119) 



4m 2 



^ 2 fc' 2 T/] 2) (a, s') + V^ 7 !^ («, a') + 8m 2 (a, a') - 8 vf (a, a') 



The following states are located on the a\ (M 2 , n)-trajectory [15]: ai(1230) with n 
ai(1960) with n = 3, ai(2270) with n = 4, and so on. 

The wave functions of the ai-states satisfy the orthogonality/normalization condition: 



' ds 



25 



p(a) fc 2 2 S 2 ^(lC'^W = a*, 



1, ai(1640) with n = 2, 



(120) 



4m 2 



2. Equation for the (13 (M 2 ,n)-trajectory 
For the a3-mesons (S = 1, £ = 3, J = 3) the Bethe-Salpeter equation reads: 



oo 

(a - M 2 ) ^A3) (a) = _ J d A p{sl) fc'6^ x 

4m 2 



(121) 



|fc 2 fc' 2 vf } s') + V^ 3) (*, + 8m 2 (s, a') - § vf } (a, a') 



Two 03-states were seen: a 3 (2030) with n = 1, a 3 (2275) with n = 2 [15]. 
The orthogonality/normalization constraint reads: 



00 

/ 

4m 2 
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P{s) ^s^a(3)^ 3 f(a)^Ha)=6 n >, n . 



(122) 



C. Equations for the (S = 1,J = L± l)-states 



We have two equations for two states with S = 1 and J = i±lforJ>0. Corresponding wave functions are 



denoted as 



(n) /J.1...HJ 



(k) + Bj* 



(n)jui. 



(fc) with j = 1,2. These wave functions are orthogonal to one another. 



Normalization and orthogonality conditions give three constraints for four mixing parameters Aj and Bj . 
Each wave function obeys two equations: 



(a - M 2 ) X^ jf} (k)Sp ^(k\ + m)^(-k 2 + m) 



(123) 
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X (^.LW^^) + fli«a(*)^' J+1 ' J) (* 8 )) = 



: W £ 5 P + ™h£(-fo + m) 



/ 



rf 3 fr' r - - - 

' ' ^^/c(s,s\(kk'))S P [^,(k[+m')F c (-k , 2 +m') 



and 



(,s - M 2 ) 41:^,^ (fc) 5p 7^1 + m) 7 ^(-fc 2 + m) 



x 1J) (k 2 ) +B J xW a {k)W>+ ^(fc 2 ) 



d 3 k' 



(2t« 



7 y c (s, a', (fcfc')) 5p 7 „'(fci + m')F c (-A£ + m A 



First, let us consider (123); in the left-hand side of (123) one has two convolutions: 

X ( ^ J0 {k)Sp k(fci + m)^(-k 2 + m)l X^ ja {k) 



2a(J)fc 2 ( J+1 ) 



2J+1 , ,,, 2 



J+l 



^SU*)^ [7a (£i + ^)7^(-fc 2 + m)l ^, a (fc) 



= 8a(J)fc 2 < J+1 ). 
The left-hand side of (124) contains also two convolutions: 

^iT.^,^^^ k(fci + m)7^(-fc 2 + m)l 4^ ja (fc) 



= 8a(J)fc 2 ( J+1 \ 



7a"(A:i+m) 7 ^(-fc 2 +m) ^.^(fc) = 



= 2a{J)k 2{J - 



2J+1 , ,,, 2 



J 



-s + 4fc / 
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The right-hand side of Eqs. (123) and (124) is determined by convolutions of the trace factor Bp/ a i (s, s', (kk 1 )), see 
Eqs. (113) and (114), with angular momentum wave functions; corresponding formulae are presented in Appendix C. 
Following them, one has for the right-hand side of (123): 



j+i 



(127) 



2J+ 1 



(V7s 7 V v (s 1 s',(kk')) + 8mm'V T (s,s',(kk'))^ + 



+As'k 2 V v {s, s', (kk')) + Ask' 2 V v (s, a') + I6^j^-k 2 k' 2 V v (a, s', {kk')) 



xP J+1 (z) + 



16mm' V s (s, s', (kk')) + 4- y/ss 1 V A (s, s', (kk')) 



and 



xVk 2 Vk l2 P J (z) + 16 ^j^T fc2fc ' 2 Vv{s,s', (fcfc'))P/-i(z) 



(128) 



= 16a(J)fc 2 (>/jfcVP) 



j-i 



s + 4 



J + 1 e 

27TT 



fc ,2 P J+1 (z)VV(s,s',(fcfc'))+ 



+ -Vss 7 V A (s,s',(fcA;'))+4mm'ys(s,s',(fcfc')) Vk 2 Vk' 2 Pj(z) + 



+ 



s' + A—^k' 2 
2J+ 1 



fc 2 Pj_l(«)W(«, «',(/!*')) ■ 



For the right-hand side of (124) one has: 



(129) 



16a(J)fc' 2 (vWifc*) 



,7-1 



s' + 4 



J+l. fc , 2 



2J+1 



fc 2 P J+1 (z)t/y( S , S ',(fcA : ')) + 



-W Va(s, s', (fcfc')) + 4mm' V s (s, s', (kk')) Vk 2 Vk' 2 Pj(z) + 



and 



+ 



s + 4 27TT fc 



Ci^w (*. ^:. 1 i,a'(fc') = 4a(j) (vpv^) 



.7-1 



(130) 



< ( 16^^fc 2 fc' 2 P /+1 (z)W( s , s ',(fcA : ')) + 
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+ 



16mm' V s (s, s' , (kk')) + V A (s, a', (kk')) 



k 2 Vk' 2 Pj(z)+ 



+ 



2J+ 1 
J 



s' (VsFVvis, s', (kk')) + 8mm' V T (s, s', (kk')fj 



+ 



+As' k 2 V v (s, s', (kk')) + As k' 2 V v (s, s', (kk'))+ 



+ 16 2jqr - I k 2 k' 2 v v (s,s',(kk')) 



Pj-i(z) 



In the right-hand sides of Eqs. (123) and (124), we expand the interaction blocks in the Legendre polynomial series 
(94) and integrate over angle variables dz/2. As a result, Eq. (123) reads: 



(s - M 2 ) 



4^-^(k 2 )A, + ^l±l s + Ak^ 4y +hJ) (k 2 )B, 



(131) 



OO 

-! 



^p(s')8(-k' 2 y-^-'-^(k' 2 )A J x 

IT 



4m' 2 



J + 1 \ l/4 T 7(J+l)/„ J 



i(J + 1)1* + J fc' 4 (5, *') + ^(J)Vss'k' 2 Vi j) (s, a') 



■"■'-<• ^) (s , s ') + £(J-i) (s' + ^j^' 2 ) ^• 7 " 1! " ■' , 



+ 



OO 

I ^p( S ')2(-fc' 2 )^Vi 1 - J+1 ' J) (fc' 2 )^x 



4m' 2 



g mm '^±l £( J + 1) V^-l4 J+1 ) ( S , S ') + 



+£(J + 1) f ^±i SS ' + 4s'fc 2 + Ask' 2 + 16^-tLk 2 k' 2 ) Vy +1) (s, s')- 



V J+l 



2J+ 1 



-16mm'f ( J) vj J) (s, s') - 4—^- £( J) Vl^vj" (s, s') + 



+16 27TT e(J_1) ^ 1)(s ' s ' ) 



The second equation, (124), reads: 



(s - M 2 ) 



B 1 



s + Ak^j ^-^\k 2 )A j+ Ak^ J ^ J \k 2 ) Bj 



(132) 



25 



DO 

■/ 



^ P (s>)2(-k' 2 y-^y-i> j \k' 2 )A j x 

7T 



4m' 2 



J 



2J+ 1 



- 1) ( ^^ss' + 4s'fc 2 + 4sk' 2 + 16 — — fc 2 fc' 2 4 J-1) ( s > 



-16mm' <£(J) fc 2 fc'Vj J) (s, a') - 4^yi £(J) vOT^V^a, s')+ 



J+l 



+16 27TT ?(J + 1)fc4fc ' 4 ^ J+1)(s ' s,) 



+ 



oo 

+ y ^ P ( S ')8(-fc' 2 ) 7 +vi i ' j+i - 7) (fc /2 )^ 



4m' 2 



e(j + 1) («' + 4 ^r fc ' 2 ) fc4 4 J+1) (*> + v^ 2 ^j J) (s, «')- 



Normalization and orthoganality conditions determined by Eq. (87) are as follows: 

oo 

/ *p(«) [A] (^■ J - 1 ' J )(fc a )) a 2a(J)(-* a )< J - 1 ) (^±l S + 4fc 2 ) + 



4m 2 



and 



+2^ J B j Vi 1 '- 7 " 1,J) (fc 2 )V'i 1 ' J+1 ' J) (fc 2 )8a(J)(-fc 2 ) ( - 7+1) 



+B 2 (^ J+1 ^)(fc 2 )) 2 2 a (J)(-fc 2 )^ (^« + 4fc 



1, j = 1, 



OO 

4m 2 



+(A 1 S 2 + A 2 i? 1 )^ 1 ' J - 1 ^ ) (fc 2 )4 1 ^ +1 - 7) (fc 2 )8 a (J)(-fc 2 )( J + 1 ) + 



+5^2 (v,„- 7+1 - 7 )(fc 2 )) 2 2 a (J)(-fc 2 ) (J+1) (jTT s + Ak '' 

Let us emphasize again: all the above equations are written for the case J > 0. 



0. 
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1. Equation for the a?, (M 2 ,n) -trajectory 



For the (S = 1, L = 1, J = 0) state, we have only one level L = J + 1; the wave function of this state obeys the 
equation: 



oo 

(s-M 2 )(s + 4k 2 )^(s) =- J ^ p{s ')2k' 2 ^\s')x 



4m' 2 



(135) 



^m 2 yfiP V™ (s, s') + £ (ss' + As' k 2 + Ask 12 + 16 k 2 k' 2 ) (s,s')- 



-16m 2 vP( S ,s r ) 

According to [15], the following states are located on the ao (M 2 , n)-trajectory: 
a (980) with n — 1, a (1520) with n = 2, a (1830) with n = 3, a (2120) with n = 3, and so on. 
The normalization reads: 



£p(«) (^W) a 2(-* a )(« + 4fc a )=l. 



4m 2 



(136) 



2. Equation for the fa (M 2 ,n) -trajectory 
The /o-states have two flavor components, nn and ss, correspondingly we have two equations for two wave functions: 

(s - M 2 ) (s + 4k 2 ) ^ } Js) cos9„ = (137) 



V^')2fc' 2 vi 1 ( 1 i ) )> „( S / )cose„ 



4m 2 



\m 2 ^V^ nn)tT (s,s')+ 



+ \ (ss' + ±-s'k 2 + Ask' 2 + 16 fc 2 fc' 2 ) V^ nn)y (s, s')- 



-^m 2 V^ nn)>s {s,s')] - J ^p.V)2l?*%$ tn W*ne n x 

Am? 



—mm s V ss v , - _s rr 

o s (ss—>nn),T 



(s, s') + \ (ss' + As'k 2 + Ask' 2 + 16 k 2 k' 2 ) x 



and 



xV, 



(i) 



(ss— mn),V 



(s, s') - 16mm s V^^ ^s, s') 



(s-M 2 ) (s + Ak 2 )^ )n (s)sine n = 



(138) 
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da' 



^p(s')2 k' 2 1>%$ )in (s') cos 0„ ^mm s Vss> V^ )tT (a, a')+ 



Am 2 



+ i (ss> + Aa'k 2 + Ask' 2 + 16 k 2 s k' 2 ) V^ sS)y (s, a')- 



— I6mm s V, 



(0) 



s (nn—>-ss),S\ 



oo 
Am 2 



W s Vs^V^ sS)T (s, a') + ^ {ss' + As'k 2 s + Ask? + 16 k 2 k' 2 ) 



x (a, a') - 16m 2 . V, w J a, a') 



2 t7(°) 



The following states are located on two /o (M 2 , n)-trajectories [15]: 

1) / (980) with n = 1, /o(1500) with n = 2, / (2005) with n = 3, /o(2240) with n = A, and so on, 

2) / (1300) with n=l, /o(1750) with n = 2, /o(2105) with rc = 3, / (2330) with n = 4, and so on. 



3. Equation for the p (M 2 ,n) -trajectory 



Two equations read: 



(s - M 2 ) 



A^\s)A + ^-s + Ak^ ^\a) Bj 



(139) 



4m 2 



3 

To 



8 + U' 4 Vf (*, *') + * V™ («,«')- 



± m 2 k' 2 V^\ SlS ')+L' + - i k' 2 \ V^(a,a') 



+ 



OO 

+ J ^ p ( s ')2k'%] n 2 ^(s')B 3 

Am 2 



m 2^ Vs^V^is, a') + ^ (Ls' + As'k 2 + Ask 12 + ^fc 2 fc' 2 ) V^(s, a')- 

±\J \ Z o 



and 
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(s M 2 ) \{3a + 4fc 2 ) ^ 1 \s)A j + Ak^^^B, 



(140) 



4m 2 



24m 2 \/ ss 7 l4 0) ( S , S ')+ ( 3s S ' + 4s' fc 2 + 4sfc' 2 + ^-k 2 k 12 ) V^ UJ ( S , S ')- 



-m 2 ^ fcW ^ (s, «') - | V^JfeWj 1 * (a, s>) + ±f fc 4 fc' 4 V^ 2) (a, a') 



+ 



OO 

+ J ^ p{s >)8k>%] n 2 ^( S >)B 3 

4m 2 



1 (V + lk'A fc 4 14 2) (5, + I V^fc 2 V™ (a,a')- 



_ m 2 4 fc 2^i) (s y )+ ^ + 4 fe2 j ^(0) (S)SO 

The following states are located on the /? (M 2 , n)-trajectories [15]: 

1) p(770) with n = 1, p(1450) with rc = 2, p(1830) with n = 3, p(2110) with n = 4, and so on, 

2) p(1700) with n = 1, p(1990) with n = 2, p(2285) with n = 3, and so on. 
Normalization and orthogonality conditions are as follows: 



OO 

J %(a) [A 2 (^„°- 1 )( s )) 2 2(3 S + 4fc 2 ) + 

4m 2 



(141) 



+2A j B j ^\a)^ 1 Ha)8k 4 + 



and 



+L> 2 (^( S )) 2 2fc 4 (^ + 4fc 2 



1, J = 1,2, 



OO 

J %(a) \a,A 2 (a)) 2 2 (3, + 4fc 2 ) + 

4m 2 



(142) 



+ (A 1 i? 2 + A 2J B 1 )^ 1) (^„ 2 ' 1) ( S )8fc 4 + 
+B 1 B 2 (^ 1 „ 2 - 1 )( S )) 2 2fc 4 Q S + 4fc 2 ) =0. 
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APPENDIX A: APPENDIX: ANGULAR-MOMENTUM OPERATORS 



Here we present the angular momentum operator X\ 11 l 2 ... liL _ ltlL (k) and briefly recall its properties; a full presen- 
tation of the angular- momentum operators can be found in [9] . 

The operator xj^l_ 2 ... f _ lL _ lML (k) is constructed by using relative momentum of mesons in the space orthogonal to 
the total momentum P: 

P V P 

kfi = k v g vil i 9vfi = Svn ~ ! 9\iv = (1> — 1j — 1j — !)■ (A-l) 

In the center-of-mass system, where P = (Po,P) = (\^s, 0), the vector k 1 - is space- like: k 1 - = (0, k). We determine the 
operator xj£j l2 ... IJ , L _ ll _ lL (k) as symmetrical and traceless. It is easy to construct it for the lowest values of L = 0, 1, 2, 3: 



= - 

(Ul/i2/i3 2 



fr 1 - fr" 1 fr" 1 - (r;- 1 A- 1 4- o 1 - i- 1 4- r, 1 - fr 1 - 



Correspondingly, the generalization of xj^}.. liL for L > 1 reads: 

X (L) = k^z {L - x) 



Z (L-1) 



2L-1 



L-D ± _ 

^i---A l i-lMi+l—ML y Mi a 



2L 



£4 

\*=i 

L 

]_ Mi— Mi-iMi+i— Mi-i^j+i- 



i,3 = l 



I 



It is seen that the operator xj£} l2 ... llL _ lftL (k) constructed in accordance with (A. 3) is symmetrical, 



X (L) 



and it works in the space orthogonal to P: 



p — n 



The angular- momentum operator x\^\. llL is traceless over any two indices: 

q 1 - X™ = 



a X (L) 

aMiMj fj,i...m...p,j ...jUl 



(A.3) 



(A.4) 
(A.5) 

(A.6) 



The tracelessness property given by (A.6) is obvious for the lowest-order operators entering (A.2), for example, 



9^ 2 X { ^ 2 = (recall that g^ 2 g^ 2 = 3). 
The convolution equality reads: 



X {L) k 1 - - (k^ 2 x {L - 1] 

Using (A. 7) we rewrite the recurrent equation (A.3) in the form: 



(A.7) 



X (L) 

H1.../J.L 



(A.6 



2L 



L 2 



1 L 

/ , Hi jiii— M.-iR+i- 



2(^) 2 
L 2 



■HL 



i=l 



■ Hj-iHj+i-.-HL 



i,j = i 
i<j 
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On the basis of this recurrent equation and taking into account the tracelessness of xj^}..^, one can write the 
normalization condition for the moment-L operator as follows: 

X^ L (k)X^.., L (k) = «(L)(^f L ,a{L) = f[ ^ = { ^-^. (A.9) 

1=1 1 ij ' 

The iteration of (A. 8) gives us the following expression for the operator xj^l. l _ iL : 



■ ■ ■ (A. 10) 



,">ii-,">ii, ■ ■ ■ ft nr T </„. * " ' ii 



2£ _ ^ (f(H(tt" /' "/'I /'/. '/' 1 /'.:" /'J " /' ; /'i. 

1 — l)(2i — 3) ^ 1 ^ 12 ' y ' J3 ' J4 ^ 5 ^ 6 ' ' ' 

One can introduce a projection operator 0£*;;££ for the partial wave with the angular momentum L. The operator 
is defined by the following relations: 

x^ L {k)o^ = x£l VL {k), 0£:::£ o«w.^ = o% ::: %. (a.h) 

For the sets of indices [ii . . . /il and v\...ul, the operator O has all the properties of the operator it is 

symmetrical and traceless, 

QHlH^—HL _ f)H2Hl—HL _ QH1IJ.2—IJ.L (QMl Ml • • -Mi _ _ A C A 1 21 

w ^lJ/ 2 ...^t w !/l!/2---^t '-'l/21'l...I'L ' w !/l!/2---^t ^1 ^1 ■ ■ -fi " V^ 1 -^^/ 

The projection operator can be consructed as a product of the operators 

x\^}.. l j, L (k)Xvi)..i/ L (k) integrated over angular variables of the momentum fc^, so we have a convolution of the (2L+1)- 
dimcnsional vectors, which provide us with irreducible representation of Lorentz group in the fe |-space. So, 

mo»\:X = j^yr I § x^ L wZLi (*), (A.i3) 

where is a normalization factor fixed below. Using the definition of the projection operator O^ 1 .'.'.'^ we have: 

k^ ■ ■ ■ ^0%;;;% = (k). (A.14) 

This equation represents the basic property of the operator: it projects any index- L operator into partial-wave operator 
with the angular momentum L. 

Multiplying the equation (A. 13) by the product X^.. f _ lL (q) ■ xif).. VL (q), we get: 

1 

Z{L)X%lvMX%l»M = (^) 2i « 2 (i) / y Pl(z), (A.15) 



-1 

that gives the normalization constant in (A. 13): 

m = ^- = / 2L - 1)!! ,. (A.i6) 

Summation in the projection operator over upper and lower indices performed in (A. 13) gives us the following reduction 
formula: 

or i i 



31 



Likewise, the summation over all indices gives us: 

0£;;;£=2L + 1, (A.18) 
that can be proven using formula (A. 13). On the basis of the equation (A. 14), one gets 

yW qui-hl-i _ y-(i) (A 19 s ) 

Generally, one can write: 

v(L) qui. ..m _ v(L) If, prA 

^fl-l.. ..fiilH+l... PL Vl—Vi V\ ...ViHi+lHL • l^.^u; 

APPENDIX B: APPENDIX: TRACES OF THE LOOP DIAGRAMS 

Here we present the traces which are used for the calculations of loop diagrams. Recall that in the spectral-integral 
representation, there is no energy conservation, s ^ s' , where P 2 — s, P' 2 — s', but all constituents are mass-on-shell: 

k-^ — — k<2 TfX , k^ — — k^ — — Tft 

The following notations are used for the quark momenta: 

K = \(h-k 2 ) v , K = ^(k[-k' 2 ) v , (B.l) 
= k v g Vjl = fc M , k'^ = k' v g v ^ = k'^. 

We follow the definition of matrices: 

75 = -H°7 1 7 2 7 3 5 <V = \ ■ 

1. Traces for the 5 = states 

For the S = states we have the following non-zero traces: 

T'p = Sp [ij 5 (k[ + m') l5 (-k' 2 + to')] = 2ia', (B.2) 
T' A = Sp [i 75 (fci + m')*7M75(-fc2 + m')] - -4m' i*, 
= 5p [i7 5 (fci + m')ia^(-k' 2 + to')] = -Aie^ a/3 P^k' , 



and 



Tp = Sp [«75(-fc2 + m)75(fci + to)] = 2zs, (B.3) 

Ta = Sp [ij 5 (-k 2 + TO)27 M 7 5 (fci + to)] = 4toP m , 

T T = Sp [ij 5 (-k 2 + m)i(T J[1I/ (fci + to)] = Aie^pPakp. 



The convolutions of the traces Ap = (TpT' p ), A A = (T A T' A ), 
A T = (TtT t ) are equal to: 



A P = -Ass', (B.4) 
A A = -16mm' (PP'), 
A T = -Z2{PP'){kk'). 
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2. Traces for the S = 1 states 



The case of the (S = 1) states, the traces are equal to: 

T' s = Sp [ 7 i(fei + to') • I • (—k' 2 + to')] = 8m'C, 

?v = ^ + m'H(-k' 2 + m')] = 2 + 4fc^) , 

= Sp [^,{k[ + m')i^ 5 (-k' 2 + to')] = 4e QW fc^, 
It - Sp [^,(k[ + m')ia^(-k' 2 + to')] = 4m'* [(^i* - <^f*] , 

and 

T s = Sp [7/3 (— &2 + to) • I • (fci + to)] = &mkp, 
T v = Sp [lp(~k 2 + to) 7m (/c 1 + to)] = 2 [g^s + 4fc /3 fc M ] , 
T A = Sp [jp(-k 2 + TO)«7 M 7 5 (fci + to)] = -lepnaipkaiPp, 
T T = Sp [yp(-k 2 + m)ia^(ki + to)] = Ami [g^P v - gp v P»\ ■ 
Corresponding convolutions B c = (T c T' c ) read: 

{B s )p' a > = Mmm'kpk'a,, 

(BvW = 4 [ss'gj, a , + 4s'kp,k a , + 4»A£,A&, + lG{kk')kp,k' a ,] , 
{B A ) p , a , = -W(PP') [k' p ,k a , - {kk')gj, a ] , 
(BrW =32(PP')mm'g^ a ,. 



(B.5) 



(B.6) 



(B.7) 



APPENDIX C: APPENDIX: CONVOLUTIONS OF THE TRACE FACTORS 

U u „± _ 



k'vSvn = an< ^ introduce: 



Here we present the convolutions of the angular-momentum factors. We work with k^ = k v g^ = k^ and k'^ = 

^ (C.l) 



(C.2) 



k 2 Vk' 2 



The convolutions for the (5 = 1) states read: 

4^ 2 .. w (fc)4^ 2 .. w (fc') = a{J) (Vk^Vk^y Pj(z). 
Analogous convolutions for (S = l)-states are written as follows: 



X 



H1H2---HJ0 



(k)X, 



(J+i) 

fJ,lfJ,2---fJ,jO! 



l(J) 



J + 



£/2 ^2 
^=^P.7,., + l( 2 ) fc/jfca + -j=B Pj J + 1 {z) k' p k' a + 

+C Pj , J+1 (z) k fj k' a + D PjJ+1 (z) k' k a + (y/V Vk^) E PjtJ+1 (z) gj a 



X 



P2 



(J-l) 

fJ.lfJ.2---fiJ,OC 



(k') = 



a{J) 1 
~l~~k Pi 



^ Ap^^z)^ + (% J+1 (z) - (2J+l)4/(*)) fc^fc^- 



+(^,, +1 «-(2J + i)£/M)M4 + 

+£> Pj , J+1 (s) *£fc Q + (vWfc*) S PjjJ+1 (z) ^ 



(C.3) 



(C.4) 
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(A PjiJ+1 {z)-(2J+l)Aj{z)) k p k 



J-2 



(B PjtJ+1 (z)-(2J+l)Aj(z)) k'pk' a + 

+ (c Pj , J+1 (z) + {2J j^' 2 Pj(z) 2(2J+l)Bj(z)j k p k' a + 
+D PjiJ+1 (z)k' () k a + (Vk^Vk^)E Pj:J+1 (z) g^ a ] , 



and 



Here, 



_ (2J + 3) 5 



a(J) (y/&y/W*} J 1 (PP r ) 



k 2 Vk^((z 2 - l)D Pj , J+1 (z) + zE PjtJ+1 (z)) gj a - 



{J + lf 



- D Pj,J+A z ) —7fW k P k ct + -7= k 'f3 k » - zk P k ' a 



+ (zDpj, j+1 (z) + E PjiJ+1 (z)) k' k o 



] (k) €ai>iV2V 3 P Vl Zv 2 Hi ■■■/J j,va ( k ) — 

= YTT^ {^^y 1 a pj.j+i( z ) t aP >kk>- 



a p.,,., + A z ) = B p.,,.,+i( z ) = 



2zPj(z)+[Jz 2 -(J + 2)]P J+1 (z) 



(1 -z 2 ) 2 
[(l-J)z 2 + (J + l)]P 7 (.) 



+ 



[(2J+l)z 2 -(2J + 3)] zP J+1 (z) 



(1-z 2 ) 



2N2 



n . , [(J + 2)z*-J]Pj(z)-2zP j+1 (z) 

^Pj,J+l{ Z ) - ^ _ z 2)2 ■ 



E Pj,J+A Z ) 



zPj(z) - P J+1 (Z) 



A,= 



(1-z 2 ) 
P J+1 (z)-zPj(z) 



1-z 2 

£ a p>kk> = eafi^P^k^k'" 



B.j = 



Pj(z) - zP J+1 (z) 



1-z 2 



We also need more complicated convolutions, namely, for the factors 
K P X ^l^WX^.. llja {k')K a where K = k,k': 



k P ^ 1 p 2 -' ) -- AI . 7 /3(^)^/ ( 'iM2-' ) -A'J«(^') k a — 

= k 2 a(J) (^Vk^Vk^Y +1 Pj+i(z), 
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^ xl ll pJ.. l _ t jp(k)X^ 1 p 2 ).. flja (k') k Q 



2J + 1 p n J p <- n 

■zPj+l(z) - — — -Pj(z) 



J+1 ' J+1 

and for the factors (^KpX ( ^l^ jp (k)z\i[~l}.^ ua (k l )K ^-. 

k P ^ 1 ^ 2 -' ) -- A ij/3(^)^iM2-' ) -MJ,a(' e ') k a = 

= fc 4 a(J) (\/fc2\/fc 7 2)' 7 ~ 1 p J _ 1 (2), 

hXli^^zlf-^Ak')^ = k 2 a(J) (y/vVw) J Pj(z), 

k 'p xX^pWztt^AV) K = «V) (vW#*) J+1 p j+1 (z), 

= fc 2 «(J) (^^) J P/(z), 

and for the factors (K/ 3 zjf i ~^^ ji p(k)Z^^jA k ') K a)- 

k P Zj llli2 }.. liJi p(tyZpiH2--HJ,a(k') k °< = 

= k 2 a(J) (V&Vu*y 1 Pj-i(z), 

k 'f) Z^nJ-./jj^ityzh^Lfij^ik') k a = 

= a{J)(y&Vk*) J ^-p-zPj^(z) - ^p-Pj(z) , 



and for the factors [K e 0VlV2V3 P Vl Zi J 2 l 1 ... llJ ^ 3 {k) e a x 1 x 2 X 3 p \ 1 



xZ (J) 

X 2 hi---ij,j 



,X 3 ( k ')K a ): 



kj3 £fiv\V2V-j,Pv\ Z V2 n 1 ...fi J >V3 (fc) £aXiX2X :i l 
kfj tfiv\V2V?,Pv\ ^v^\-l—\ij,vz (k) e aXiX 2 X 3 P\ 1 Z 
kp (■jiv\V2V 3 Pv\ ^ii^i\—\ij,V3 (k) e aXiX 2 X 3 P\ 1 Z 



(2J + 3) 

-a 



AJ) 

'V2V\ 

2 , , , sj+l 



(J) 

X2H1- 


■HJ,X 3 ( k "> ka 


= 0, 


(J) 

X 2 fJ-i- 


■■^j,X 3 ( k ') k 'a 


= 0, 


(J) 

X 2 fJ-i- 


■■Hj,X 3 ( k ') k 'a 


= 0, 


(J) 

X2H1- 


■nj,x 3 ( k '} ka 





(J)(vWfc«) (PP')[zPj(z)-P J+1 (z)}, 



(J +1)3 

9(iu t l3v\V2V 3 Pv 1 Z V2 y L1 ...y LuV?i {}i) ^aXiX2X 3 Px 1 ^X2Hi--HJ,X 3 ( k ) = 

-a(J) (VFVfc^) J (PP')Pj(z). 



J(2J + 3) 2 / — • — W 



(J +1)3 
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APPENDIX D: APPENDIX: THE BETHE-SALPETER EQUATIONS FOR THE w, <f>, A 2 AND F 2 

TRAJECTORIES 



Here we present the Bethe-Salpeter equations for the oj, <f>, a 2 , and / 2 trajectories. Though the explicit form of 
these equations is rather cumbersome, the investigation of these trajectories is informative for the reconstruction of 
quark-antiquark forces. 



1. Equations for the a 2 (M 2 , n)-trajectories 



The following states are located on two (n, M 2 )-trajectories for the a 2 states M < 2400 MeV [15]: 

1) a 2 (1320) with n = 1, a 2 (1660) with n = 2, a 2 (1950) with n = 3, a 2 (2255) with n = 4; 

2) a 2 -trajectory: a 2 (2030) with n = 1, a 2 (2310) with n = 2. 
Correspondingly, we have two coupled equations for two wave functions: 



(s - M 2 ) 



(D.l) 



and 



00 

■/ 

4m 2 



da' 

^ p(s 08(-fc' 2 )^ 2 )( S ')^x 

7T 



- V k' 2 vP(s, s') + Us'+ °-k< 2 ) V^(s, s') 



00 

-/ 

4m 2 



da 1 

^p(s')2k'^ rl 2 \s')B j x 



m^V^V^is, a') + A ( + 4s'fc 2 + 4sfc' 2 + ^fc 2 fc' 2 ) x 





21 



14 V 3 



x (s, s>) - m 2 I > (5, - 1 V^Vf *') + I V™ (s,s>) 



(s - M 2 ) 



^.s + 4fc 2 ^{s)A 3 +Ak^\s)B 3 



(D.2) 



DC- 

■/ 

Am 2 



d J- p{s >)2{~k> 2 )^X s i)A 3 * 



_ m 2 ^ fc 2 fc' 2 ^ 2 ) ( S , «') - £ vOTfc'Vj 2) ( a , a ') + ^ A: 4 fc' 4 v4 3) (a, *') 



3G 



ds' 



Am 2 



14 



12. 



77 I s ' + ^ fc4 V V (*> *') + TTT Vi7fc 2 (a, S ') 



10 



Normalization and orthogonality conditions read: 

ds 



J ^p(s) [A 2 (rt\¥\s)) 2 2a(2)(-fc 2 ) (jj, + 4fc 2 ) + 

4m 2 



and 



+2A j B j ^\s)^\s)8a(2)(-k 6 )+ 



+^(V& 8 „ lJ) 00) 2a (2)(-fc 6 )f| s + 4fc 



= 1 



ds 



P(s) 



Am 2 



MAiB 2+ A 2 B 1 )^\ S )^\ S )%a{2)(-^) + 



+B lB2 (^ 3 „ 2) ( S )) 2 2a(2)(-fc 6 ) (h + Ak 



= 0. 



2. Equations for the uj and cj> (M 2 , n)-trajectories 

We have four trajectories in this sector with the following states located on the (n, M 2 )-trajectories [15] 

1) .S-wave dominant states: u(780) with n = 1, w(1420) with n = 2, w(1800) with n = 3, w(2150) with n = 

2) D-wave dominant states: w(1640) with n = 1, w(1920) with n = 2, w(2295) with n = 3; 

3) S'-wave dominant states: 0(1020) with n = 1, 0(1660) with n = 2, 0(1950) with n = 3; 

4) £>-wave dominant states: 0(1700) with n = 1. 

Correspondingly, we have four coupled equations for four wave functions. The first one reads as follows: 



(s - M 2 ) 



oo 
4m 2 
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m 2 k' 2 V, W 



, 4 /.'2^ t7(0) 



+ 



oc 

+ / ^P( s ')2fc'Vi 1 ; ( 2 i ) ) ,„( s ')^cose;x 

4m 2 



i,r^v'«'^ _^( S , S ') + JL (| ss ' + 4s 'P + 4 S fc' 2 + yfc 2 fc' 2 J 



,18 



5 V (nn — mn),T*> ' ' 10 \ 2 



' (nn—>nn),V ^ 



< 1 ; - , „(,, S ') - m 2 ^ V^ nn)iS (s, s') - I ^V^ nn)iA (s, s> ) + 



o (nn— *nn),V v ' ' 



+ 



/ v^ (s ' )8 ^(S» (s,) ^ sin0 " x 



4m 2 



+ 



OO 

-J ^ Ps ( s 02fc;v^„( s ')^sine;x 

4m 2 



mm »y^^„ s)iT (s, s ') + ^ Q ss ' + Vfc 2 + 4s fcf + |fc 2 fcf) / 



xV, 



(2) 



(s, - mm^^^,^, s') - H V^V^ n)iA {s, s')+ 



16 



+ 3 ^(is ) -^nn),y( S ' S ') 



the second one, 



(s - M 2 ) 



sine « + (r + 4fc2 ) <S ) ,»w x 



(D.6) 



OO 

cosine;] = J ^ P ( s ')S^ )n { S ')A jC osQ n , 

4m 2 



3 / 8. 
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the third one, 



\mm s k' 2 v£L„-),s(*. <0 + + f *") ^°L ss -) ; v 



(*,*') 



+ 



OO 

h | ^ p(s ' ) 2fc'Vi 1 fi ) ) ,„( s ')^cose;x 

4m 2 



18 
"5 



32.0 
3 



( ( „ 2 L ss -),y( s ' s ') - mm «y ^ ( r!L« s -),s( s ' s ') - 3 ^!!L 8S -),a( s . s ')+ 



+ f s s >) 



4m? 



1Q \ O SIS {SS^SS),V V ' / O 8 (s.S—>-SS),A^ ' / 



+ 



OO 

4m? 



,18 



„a ,3 /3_ , 



10 V 2 



32 



16 



3 ^(ss— s g),y( s ' s ') 



( S - M 2 ) I (3s + 4fc 2 ) <Jii lB («)^ cosQ ™ + 4fc4 <s;,„(«)x 



4./,(1.2,l) 



(D.7) 



OO 

xB jCO se;]= I ^p( s ')2<fi ) ) ,„( s ')^cose„) 

4m 2 



1 .s .s' ) + ( 3ss' + As'k 2 + Ask' 2 + y fc 2 fc' 2 ) , 



. r i; ; WA ( S;S ')- m 2_^fc 2 fc' 2 K (1) 
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« {s,s> ) + l 4kwm nYV {s,s>) 



(nn—>nn),A\ ' 5 (ran— >rara),V^ 



ds' 



4m 2 



X Bj COS 0^ 



1 ( s > + 8 fc ' 2> ) fc 4 vf ) _. _,(«, a ') + 1 V^fc 2 , («, a')- 

\Q \ 3 / (ran— »ran),y\ ' 3 (ran— »rara),AV"' / 



2 \ t7(o) 



+ 



+ 



I ^ Ps ( S ')2<' ( °i ) ) ,„( S ')^sine n x 



4m 2 



24 mm, V ss' V, 



(o) 



* {ss^nn),T 



(.s. i + ( :!,,' + ls'/, J + I,// 2 + yfc 2 fc^ ) 



16 



2i./2t>(1) 



(ss^ran),V 



3 s (ss^nra),S>- 



M .1^1.4, /4 T 7(2) 



4m 2 



x sin 9^ 



-mm s -e v£> s (8, s')+[s+ ifc 2 ) V£> v (8, s>) 



3™ K (ss— >rara),SV ' ° / ^ \ 3 / ( sS ^ ra ")' y 



2\ f>(0) 



and the fourth equation: 



(s - M 2 ) I (3s + 4fc s 2 ) sin6„ + 4*>™„( S )x 



4 ,,(1,2,1) 



(D.8) 



OO 

4m 2 



•-'-'.<""' - n " : *'V[°l^ sS)tT (s,s')+ (3 SS ' + 4 S 'fc 2 +4 S fc' 2 + yfc 2 fc' 2 ) < 



16 



.2j./2f>(l) 



(rara— >ss), V 



/i.2,/2 T 7(l) 



c s (rara— »ss),v \ ' ' 



+ / v^' )8fc,4 <S«).« (s,) ^ cose « x 



4rra 2 



40 



A f s ' + !tisA /.I r-'-'i 



i I./7^,2 T 7(1) 

3 



s- + -k rz ) kt v>'L As, s') + - Vss'k'i vyi As, s'y 



(nrt— >ss),V 



+ 



oc 

+ J ^p s (s')2^l n (s')A jS me n x 

Ami 



24m^^ ss - );T ( S , s') + (sss' + 4s'fc s 2 + Ask? + ^fc s 2 fcf ) * 



xV, 



&*),v(* - ™2 y tftf^^k s') - | ^kXXls^, s')+ 



■ —k 4 k' 4 V (2) (s s'l 



4m 2 



-m 2 -fc 2 V/^ „(*, s') + ( a + -kl ) V'"l «') 



2 t7(0) 



3. Equations for the /2 (n, M 2 )-trajectories 



The /2 mesons lay on the following four trajectories in the (n,M 2 ) plane [15]: 

1) dominantly P-wave states: / 2 (1285) with n = 1, / 2 (1640) at n = 2, / 2 (1950) with n = 3, / 2 (2210) with n = 4; 

2) dominantly P-wave states: / 2 (1525) with n = 1, / 2 (1790) with n = 2; 

3) dominantly F-wave states: / 2 (2020) with n = 1, / 2 (2290) with n = 2; 

4) dominantly F-wave state: / 2 (2200) with n = 1. 
We have four equations for these four sets of states: 
the first one, 



(s - M 2 ) 



(D.9) 



x P., cos 



]= / V P( * ,)8( " fc ^A:S),n(O^COSe n X 



4m 2 



A (a + H k A k» V^ nn)y (s, s') + ^fc' 2 V^ nn)tA (s, s')- 



1 



y. ° \ T 7(i) 
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4m 2 



100 
~21 



m 2 ^^^L BS) , r (v')+ 



5 /5 



- m 2fly( 2 ) _ ( a s ') _ _ y^/V^ -n .(s, s') + — V/ 



24 

^ v (nn—>nn),S\ 



(*,*') 



+ 



oo 
4m? 



3 



{8,8') 



OO 

- J ^ Ps (8')2kf^ ) j8')B j sme' n x 

4m 2 



the second one, 



15 (sS^rm),yV 6 ' 6 J 



(s - M 2 ) 



(D.10) 



OO 

xB jS mG' n }= J ^ p ( s ')8(-fc' 2 )^i 2 4„mcose„x 

4m 2 



5/12 
14 V 5 



^ ( 8 + ^ ) V$l^ v (8, 8') + ^k' 2 VgL^ tA (8, 8') 



■%mm«k' 2 V^l ^ Js,s') + Us'+ h' 2 ] V ( {1) 



5™"'"' (nn^ss).S^ 7 3 I 5 / (nn^ss),V 
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4m 2 



100 



'^L„-),r(*.*')+ 



5 (h 



+ i4 1 i ss ' + 4s ' fc * + Ask ' 2 + A i k ' k ' 2 ) ^("-)^ s < «' 1 



24 7(2) 



( S , S ')-5^<L SS -),^ 



+ 



oo 
4m? 



-> 2 *? itf5L«. s ('. + £ ( + ^ 2 ) vA 1 '.^ 



'^1 (.>±*V* \ t7(!) 



oo 
4m? 



2 100 ^-7.7(3) 



(S, S') + A Q SS ' + 4 S 'fc 2 + 4 S ^ 2 + |^ 2 ) 



15 (ss^ss),yV*; 6 J 



the third one, 



(a - M 2 ) 



(D.ll) 



OO 

xSicose;]= | ^ p ( s ')2(-fc' 2 )^ ; ; n 2 ) )n ( s ')^cose„x 

4m 2 



3 W 



32, 
~5~ 



xK^ -, v (a,s')-m 2 — fc 2 */ 2 ^ Js, s') - - V7?k 2 k' 2 V ( {2 l , 



Hi) 
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+— k^V^ (i sM 



ds' 



^p( s ')8fc' 6 ^;g );n (s')^cose;x 



4m 2 



12, 



+ 



oo 

+ J ^ Ps ( s ')2(-^) s ^;S ) ) ,„( s ')^sine„ 

4m? 



20 / (i) 

mm s — Vss' V {sg ^ n - ) T {s, s ) + 



24 



„2./2 x/ (2) , /x _ 9 j— 7,2^/2^(2) 



-™». y W^L^sO*. s ') - ^ ^ fc2 s')+ 



(ss— >nn),A ^ 



-—k 4 k' 4 V {3) (s s'l 



oo 

- 1 ^ Ps ( s ')8^S,„(«')^sine; ; 

4m 2 



A ( s ' + l^fef ] A: 4 K (3) v (s, s') + — V7?k 2 V™ a')- 



and the fourth equation, 

( a _ m 2 ) [g S + 4k^j ^iMAi sine„ + 4fc>g;g ;n ( S )x 



(D.12) 



OO 

cosine;] = I ^( s ')2(-fc' 2 )^;;„ 2 4„( s ')^co S e n x 



4m 2 



mm, 



3 V 2 



24 



.2l.'2f>(2) / /n _ 9 ,/-T7 i ,2w2 T 7(2) 



^ P ( s ')8fc'Vg;S );n ( s ')^cose; ; 



4rn 2 
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, r -//- I /, ' 1 ' " „(*, S ') + — \/i?fc s 2 AS, s') 

' s (nn—>ss),V\ ' ' 1Q s (nn— >ss),A^ ' / 



ds' 



im1 



20 



+1 Q SS ' + As'kl + 4< + §^ 2 ) V { ^ )y (s, s')- 



_ m 2 ?1 k 2 s k'M 2 } «.(«, «') - - VZPkZtfV™ , (a, «')+ 

8 c 8 8 (SS—>SS),S y 1 ' C s s (ss— >Ss),.A v - ' / 



24 



- /rV 4 T/ (3) (<t <t'} 

7 K S K S V ( s g^ s s),vl s ' s ; 



^0^,(1,3,2) ^D fSin e' x 



-Ps(sl8k:^ )tn (s') Bji 



4m 2 , 



o 
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FIG. 1. Nonhomogeneous Bethe-Salpeter equation for the scattering amplitude; dashed block is the interaction kernel. 




(d) 

FIG. 2. a) Presentation of the scattering amplitude as a set of the ladder diagrams with the f-channel meson exchange 
interaction; b,c) cuttings of the ladder diagrams; d) meson production processes which are determined by ladder diagrams. 
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(b) 




FIG. 3. a) Scattering amplitude near the pole related to the bound state; b) Bethe-Salpeter equation for the bound state 
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FIG. 4. The partial-wave amplitude singularities in the complex-s plane. 
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FIG. 5. Scattering amplitude in the dispersion-relation approach as a set of loop diagrams with separable vertices; there 
no energy conservation in the intermediate states, s =fi s' ^ s" , and so on. 
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